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Abstract. 

Carbon-based nanostructures and graphene, in particular, evoke a lot of interest 
as new promising materials for nanoelectronics and spintronics. One of the most 
l2 \ important issue in this context is the impact of external electrodes on electronic 

r^ ' properties of graphene nanoribbons (GNR). The present theoretical method is based on 

lyj . the tight-binding model and a modified recursive procedure for Green's functions. The 

D I results show that within the ballistic transport regime, the so called end-contacted 

geometry (of minimal GNR/electrode interface area), is usually more advantageous 
for practical applications than its side-contacted counterpart (with a larger coverage 
area), as far as the electrical conductivity is concerned. As regards the giant 
magnetoresistance coefficient, however, the situation is exactly opposite, since spin- 
'■^ ' splitting effects are more pronounced in the lower conductive side-contacted setups. 

o 
o 



PACS numbers: 81.05.ue, 75.47.De, 73.23.Ad 
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1. Introduction 

Different contact configurations affect the conductivity, and thereby also giant 
magnetoresistance (GMR) in the case of ferromagnetic electrodes. Here we study 
two most common geometries, viz. the end-contacted (EC) geometry and the side- 
contacted (SC) one. In the former the interface area is minimal, whereas in the 
latter a finite coverage area exists. Theoretically, it is easier to work with EC setups 
rather than, experimentally more relevant, SC ones. The former approach might 
be justified when a nanostructure/electrode coupling is strong (good transparency of 
interfaces, weak chemical bonding). There are experimental evidences supporting this 
view. In particular, Mann et al. [1] were able to control a length of CNT (carbon 
nanotube) under a metal contact and showed that ballistic transport in Pd/CNT 
occurred basically at the edges. Moreover, in similar systems, no observable difference 
in I-V characteristics between side-contacts and end-contacts was found [2]. In general, 
however the contact configuration does matter; as shown in [3] for Pd and Ti contacted 
carbon nanostructures. Likewise, ab initio studies in [1] show that in the case of ultra- 
small graphene nanoribbons, the EC configuration outperforms the SC one, since a large 
interface area may worsen the contact conductance. 

As concerns ferromagnetic electrodes, it has been shown that in order for a 
transition metal/graphene junction to have good spin transport properties it is necessary 
that no strong chemical bonds (engaging graphene's tt electrons) are formed [S], E]. 
Although in the case of typical ferromagnetic electrodes this requirement is hardly 
fulfilled [7], adsorption of graphene (e.g. on Ni or Co) can be effectively weakened by 
introducing a protective interfacial monolayer of Cu [S]. Noteworthy, high-performance 
ferromagnetic-metal/graphitic-nanostructure contacts have already been reported for 
Nii_a.Pda; with X ~ 0.3 -r- 0.5 (see [S] and references therein). 

Although the EC geometry is much more frequently used by theoreticians 
just for reasons of computational convenience ([Sj-[I1]), the issue of side-contacted 
nanostructures has already also been raised, see [21 IH [131 HHj (paramagnetic electrodes) 
and [HI [71 [16] (ferromagnetic electrodes). In this study we treat both the geometries 
on the equal footing, and show how the corresponding transport characteristics 
(conductivity, shot noise Fano factor and giant magnetoresistance) compare with each 
other. 
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Side-contacted geometry (SC) 
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End-contacted geometry (EC) 

Figure 1. (Color online) Schematics of the side-contacted and end-contacted 
setups. S and D denote the source and drain electrodes. The GNR is drawn 
in black, and W, L, Lq and Rq stand for its width, suspended length and 
electrode-supported lengths, respectively. 

2. Method and Modelling 

We study graphene nanoribbons with infinite S-diinensional metallic electrodes. Details 
on modelling and the computation method may be found in [HI [T7]. However, for 
completeness, the most important definitions are listed below. 
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GMR = 100 1 
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Here g is the conductivity of the GNR, having width W and length L, for the parallel 
or antiparallel alignments of electrode magnetizations {ss' =tt) ti? respectively). Tss',s", 
in turn, stands for the transmission matrix of s"-spin electrons in the ss' configuration. 
Finally, GMR is a giant magnetoresistance coefficient. 
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The modified recursion procedure for SC setups is carried out stepwise. Tlie 
vertical contribution to tlie transmission matrix is related witli tlie self-energy E;^ = 
Tio,is Qs Tis,io, where I=L (R) is the left (right) electrode index, IS (10) enumerates 
electrode (graphene) interface atoms and gs is the surface Green's function matrix of the 
electrodes. The starting local Green's functions are defined as g^ = {E — Dlq — Eg)~^, 
g§_^.l = {E — Drq — Sf )^\ where LO (RO) are the electrode-supported parts of the 
GNR (with the matrices of hoping integrals D), whose suspended part extends from 
the periodicity unit number 1 up to N. Figure 1 schematically shows the systems under 
consideration, and defines the end-contacted and side-contacted geometries. 

The calculations have been carried out using a single-band (vr orbital) tight- 
binding Hamiltonian, with the GNRs strongly coupled to infinite ferromagnetic external 
electrodes (assumed to be 50% spin-polarized), and the adopted approach is based 
on the Landauer type formalism. As mentioned in Seel quasi-ballistic transport in 
carbon nanostructures has been reported in numerous papers (e.g. [in],[IS|-|2I])- Here 
we ensure a good transparency of the GNR/ contact interface by defining the hopping 
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Figure 2. (Color online) The metallic zigzag GNRs and the semiconducting 
armchair GNRs of comparable sizes (W ~ 8nm, L ~ 22nm, and A~l/3), with 
the indicated electrode-supported lengths Lq: conductivities (main panels), 
GMR coefficients (left Insets), and shot noise Fano factors (right Insets). The 
parallel (antiparallel) electrode alignment is marked with Yt iti)- 
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integral across the interface as a geometrical mean of the respective hopping integrals 
of the GNR and the electrode (c/. P El). 



3. Discussion of the Results 

In the SC case the transport properties depend on the interface area. For small area 
interfaces the results are similar to those corresponding to the EC geometry, otherwise 
the results converge with the increasing area. In what follows, the SC results, which will 
be presented, correspond to the GNR/electrode coverages large enough to guarantee the 
convergence. The representative results are presented in Figs. 2-4, as a function of the 
GNR Fermi energy which is experimentally controllable by gate voltage. The aim is to 
elucidate the role of the interface character (EC vs. SC), the aspect ratio (A=W/L) 
and the current flow direction (zigzag vs. armchair) on spin-dependent transport in 
2-terminal setups with ferromagnetic contacts. The convention adopted here is the 
following: solid lines apply to the EC geometry, dashed lines - to the SC geometry, 
whereas the cyan (brighter) curves and the black ones correspond to the parallel and 
antiparallel alignments of the ferromagnetic electrodes, respectively. The figures 2-3 
correspond to the metallic zigzag GNRs and the semiconducting armchair GNRs of 
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Figure 3. (Color online) As Fig. [2] but for the short systems with the aspect 
ratio A~3.5. 
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Figure 4. (Color online) Spin-dependent conductivities, GMRs and Fano 
factors for metallic armchair GNRs in the end- and side- contacted geometry. 



comparable sizes. In shorthand notation the systems in question are: zzl9-90 and 
ac35-52 (Fig. 2), and zzl9-10 and ac35-6 (Fig. 3), where the first number denotes the 
GNR's width and the second one (after the hyphen) refers to the suspended length, 
all these numbers are expressed in periodicity units appropriate for the given direction. 
Remarkably, for a sufficiently long coverage length (Lq) the long zigzag GNR reveals 
a saturation effect as concerns the spin- dependent conductivity (at conductance= A- 
conductivity ~ 2e^//i), resulting in gradually vanishing of both GMR and F. These 
features may be understood, as due to development of a low-energy perfectly conducting 
channel with no backscattering counterpart (c/. a paramagnetic electrode case of |17j). 
Interestingly the results for the EC and SC configuration do not differ much in the 
zigzag case; in other cases the respective differences are larger, unless the electrode/GNR 
overlap is substantially smaller than the indicated SC values of Lq in the figures. It has 
been recently shown by ab initio calculations that Ni/GNR/Ni setups, with very small 
overlap lengths (up to 8 A), have almost identical contact resistivities which moreover are 
independent of the current orientation [16j. It should be stressed in this context that the 
present phase-coherent transport theory overestimates interference processes, which in 
reality may be suppressed by diffusive and dephasing mechanisms. Anyway, within this 
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approach the conductivity does depend on the current propagation direction unhke [16] 
but in accord with [221 123] • Most probably this discrepancy conies from different aspect 
ratios of the studied systems. Noteworthy, the present theory correctly reproduces 
the so-called universal minimum conductivity value predicted theoretically [21] and 
confirmed experimentally [201 ES], since the low-energy conductivities corresponding 
to the high aspect ratio EC systems are close to 4e^/(/i7r) (see Figs. 3-4). 

Figure 4 shows the results for nominally metallic armchair GNR ac34-52 and ac34-6 
of different length. The former (longer) reveals pronounced Fabry-Perot oscillations, 
resulting from the linear energy spectrum in this case. Both the systems behave 
similarly to the other armchair systems discussed so far, in particular the conductivities 
corresponding to the ES geometry are clearly higher than those referring to the SC 
geometry. At the same time, the EC Fano factors are, on average, lower than the SC 
ones. So, the EC geometry seems to be optimal for nanoelectronic applications. However 
the accompanying spin effects - spin splitting of the conductivities and the GMR - are 
most often the strongest in the side contacted systems. 

The difference in conductances between the EC and SC configurations provides 
information about a contact resistance. If both the conductances are roughly equal 
to each other than the contact is characterized just by the channel width, otherwise 
the areal component comes into play [7|. It turns out that the contact resistances 
may be vanishingly small, especially at gate-controllable resonances (see Figs 2 and 
4). Incidentally, very small contact resistances have been experimentally found, as 
well, in nearly ballistic graphene/metal junctions [15]. In the case of more complicated 
electrodes than those studied here, extra contributions to the contact resistance can 
be important (e.g. spreading resistance [261 EZ] )• It is so, when an STM (STS) tip 
is used as an electrode and/or there is a constriction in the system. In the context of 
spintronics related studies of graphene, it should be also mentioned that apart from 
the GMR also the spin-torque-transfer (STT) problem is crucial for magnetic memory 
applications. These two akin phenomena are due to a long spin-relaxation length in 
carbon nanostructures and, in principle, the present method has a potential to describe 
the STT, too; obviously when one of the electrodes is a non-massive free layer (see 
[221 [231 EH] for SST in EC-type setups). 

4. Summary 

In conclusion, the present quasi-ballistic theory shows that, spin-dependent transport 
through GNRs, depends on their aspect ratio, the GNR/electrode interface, and the 
current direction. It turns out that the conductivity of the semiconducting GNRs, 
as well as that of the nominally metallic armchair GNRs, is clearly higher in the EC 
geometry than in the SC geometry. However in the case of long zigzag GNRs, which 
are always metallic, both the geometries yield comparable conductivities. As concerns 
the GMR coefficients and the Fano factors, they are always anticorrelated with the 
conductivities, and typically the SC GMR effect is higher than the EC one. Moreover, 
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practically in all the cases the range of the appearance of the so-called inverse (negative) 
GMR is strongly reduced in the SC geometry. 

The new insight provided by this study is that high aspect ratio (wide and short) 
GNRs appear good candidates for interconnects in miniaturized electronic devices due 
to their preferential EC (edge) type conductance. As concerns miniaturization in 
spintronics, the short EC armchair GNRs, with GMR~ 20% seem also quite attractive. 
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